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Abstract 

I introduce two new boundary layers to the optimal domain-wall 
fermion action, and formulate vector gauge theories (e.g., QCD) with 
the new action. The generating functional for n-point Green's function 
of quark fields is derived. The bilinear and higher-order quark observ- 
ables are shown to be as-invariant, to possess optimal chiral symmetry 
for any finite N s , as well as to preserve continuum discrete symmetries 
at finite lattice spacing. With 2D abelian background gauge fields, I 
compute the index, the axial anomaly, and the fermion determinant 
(ratio) of the effective lattice Dirac operator, for the optimal DWF and 
the conventional DWF respectively. All quantities of the optimal DWF 
are in good agreement with the exact solution, but the fermion determi- 
nant (ratio) of the conventional DWF disagrees with the exact solution, 
for topologically nontrivial sectors. 
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1 Introduction 



The basic idea of domain-wall fermions (DWF) E| is to use an infinite set 
of coupled Dirac fermion fields [ ip s (x),s G (—00, 00) ] with masses behaving 
like a step function m(s) = m8(s) such that Weyl fermion states can arise as 
zeromodes bound to the mass defect at s — 0. However, if one uses a compact 
set of masses, then the boundary conditions of the mass (step) function must 
lead to the occurrence of both left-handed and right-handed chiral fermion 
fields, i.e., a vector-like theory. For lattice QCD with DWF f3J, in practice, 
one can only use a finite number (N s ) of lattice Dirac fermion fields to set up 
the domain wall, thus the chiral symmetry of the quark fields (in the massless 
limit) is broken. Therefore, a relevant question is how to construct a domain- 
wall fermion action such that the effective 4D lattice Dirac operator has the 
optimal chiral symmetry for any finite N s . Nevertheless, the chiral symmetry 
is not the only concern of lattice QCD with domain-wall fermions. Another 
crucial question is whether the effective 4D lattice Dirac operator depends on 
the lattice spacing in the fifth dimension. If the dependence on is strong, 
it is unclear what is the proper value of 05 for measuring physical observables, 
assuming N s and the lattice spacing a are fixed. 

Recently, a lattice domain-wall fermion action has been constructed [5] 
such that the effective 4D lattice Dirac operator preserves the chiral symmetry 
optimally for any finite N s , and is invariant for any 05. Further, it has been 
asserted that the effective 4D lattice Dirac operator for the internal fermion 
loops is exponentially-local for sufficiently smooth gauge backgrounds jB]. 

Explicitly, the action of optimal lattice domain-wall fermions reads 1 

N s 

Af = J2ip(x,s){(u s a 5 D w (x,x f ) + 5 XtX/ )5 S:S/ 

s,s'=l x,x' 

-(5 XtX > - u s a 5 D w (x, x'))(P_5 s > >s+1 + Pj r 5 s >^ 1 )}ip{x l ', s') 

s,s'=l x,x' ~~ 

with boundary conditions 

P+tp(x, 0) = -rm q P + ij(x, N.) , r = - — (2) 
P^(x,N s + l) = -rm q P-if)(x, 1) , (3) 

where m q is the bare quark mass, D w is the standard 4D Wilson-Dirac operator 
plus a negative parameter, 

4 

D w = Yl^ + W-mo, m e (0,2) (4) 

M 

this paper, we suppress the lattice spacing a, as well as the Dirac and color indices, 
which can be restored easily. 
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t^x.x') = - [u iM (x)6 x > >x+lx - Ul(x')8 x > iX -^ , (5) 
W{x,x) = ^-\28 x ^ -U^(x)5 x ^ x+ ^-Ul(x')5 x ^ x -^ , (6) 



and the weights 



u s a 5 = — ^ — \Jl - k /2 sii 2 (v s ; k'), s = 1, • • • , N s . (7) 
Here sn(w s ; k') is the Jacobian elliptic function with argument v s f|12|) and 



modulus k' = Jl — ^min/^maxy where \ m ax and A mi „ are the maximum and 
the minimum of the eigenvalues of l-H^I, 

H w = l5 D w . (8) 

The weights in (JJJ) are obtained from the roots (u s = (u s a 5 )~ 2 , s = 1, • • • , N s ) 
of the equation 



Sz(u) 



1 - y/uR { z' n) (u) = , N s = 2n + 1 
1- ^uR { z' 1,n \u) = , N s = 2n 



such that the operator S op t (122)1 in the effective 4D lattice Dirac operator is 
equal to the optimal rational approximation of sgn(H w )=H w (H 2 )~ 1 / 2 , 



q o rn ) ) H ™ R z' n \ H l), N s = 2n+1, 
bopt ~ bopt{Hw) ~ 1 H w Rt^ n) {Hl\ N s = 2n. 



where Rz(H^) is the Zolotarev optimal rational polynomial [7||H] for the inverse 
square root of H^, 

Amin i =1 1 + H w l C2/-1 

and 

p (n-l,n)/ tt2\ _ ^0 II"=i 1 (l + h^/c'21) , v 

Amin lh=il i + n wl c 2l-\) 

and the coefficients <i , <ig, Q and are expressed in terms of elliptic functions 
jH] with arguments depending only on N s and A^ aa ./A^ n . 
The argument t> s in (JJJ) is 
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where 



I 2 



sn 2 (i 2 ^;*/) 



X' is the complete elliptic integral of the first kind with modulus k', and G 
is the elliptic theta function. From (J2J), it is clear that \~} ax < a 5 cu s < A~* n , 
since sn 2 (; ) < 1. 

Obviously, the optimal domain-wall fermion action ((TJ) is invariant for any 
a 5 , since its dependence on a 5 only through the product u s a 5 (JTj) which only 
depends on N s , X min and tfnax/^min- Therefore, a 5 is a redundant parameter 
in the optimal domain-wall fermion action ((TJ). 

Now the effective 4D lattice Dirac operator of the optimal domain-wall 
fermion action can be extracted from its fermion determinant 

detP(m g ) = J [dip)[dip)e- A f [ ^ ] , (15) 

where T>(m q ) is the 5-dimensional optimal DWF operator defined in ((TJ). The 
Grassman integral on the r.h.s. of ()15|) can be evaluated by transforming 
{4>s3s] to {f]s,Vs} according to 

01 = P-7]i + P + T] 2 , 0i = fjiiuj^Hw - 1) _1 75 , (16) 
ip s = P-7] s + P+i] s+1 , ip s = fj s (uj s a 5 H w - 1) -1 75 , s = 2, ■ ■ • , N s - 1 (17) 
i>N s = P-Vn s + P+Vi > ^iv s = f] Ns (u Ns a 5 H w - 1) _1 75 , (18) 

and then integrating {T] s ,fj s } successively from s = N s to s = 1. The result 
(will be derived in Section 2) is 



det£>(m g ) = if det[r(D c + m,)] (19) 



where 



K = l[det[ l5 {u s a 5 H w -l)}det[-P + + l[T- 1 P_) , (20) 



8=1 \ S = l 



r£ c = i±^, (21) 

i — n N i t s 

So Pt = 1 + n g T ; > ( 22 ) 

T s = \^§L , s = l,-,N.. (23) 
l + uj s a^n w 
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Even though D c is chirally symmetric in the limit N s — > oo, one still 
cannot identify (D c + m q ) as the effective 4D lattice Dirac operator for the 
internal quark fields, since D c is nonlocal and det(D c + m q ) does not reproduce 
det[7 At ((9 /J + iA^) + m q ] in the continuum limit (a — > 0). (It is understood that 
the fermion determinant is normalized by its free fermion limit.) Also, one 
cannot neglect the factor K which can be written as the determinant of a 4D 
lattice fermion operator. 

The resolution is to introduce a set of pseudofermions (Pauli-Villars) fields 
{(f), 4>} similar to the domain-wall fermions but obeying Bose statistics. Ex- 
plicitly, the action of the pseudofermion fields reads 

N s 

A p f = s){(uj s a 5 D w (x,x') + 5 XtX >)5 StS > 

s,s'=l x,x' 

-(S x y - u s a 5 D w (x,x'))(P-5 s > >s+1 + P + (5 s / jS _i)}0(x', s') (24) 
with boundary conditions 

P+0(x, 0) = -P+<j>{x, N s ), P_0(x, N s + 1) = -P.<f>{x, 1) , 
where m q = 2mo has been imposed. Then, we have 

det[D(m q )V-\2m )] = /[#][#] [#][#] e -^^-^^ 

= det[r(P c + m 9 )(l + rP c )- 1 ] . 

Thus the effective 4D lattice Dirac operator for the internal quark fields (dressed 
with pseudofermions) can be identified as 

D{m q ) = (D c + m q )(l + rD c y l = m q + (m - m q /2) (1 + j 5 H w R z {H w )) (25) 

where Rz{H^) is the Zolotarev rational polynomial defined in The 
exponential locality of D(m q ) (for any m q and N s ) has been asserted for suffi- 
ciently smooth gauge backgrounds jHj. 

In the limit N s — > oo, S opt becomes the sign function of H w 

lim S opt = 



and D c (|21j) is exactly chirally symmetric 



lim (P c75 + 75 P C ) = , (26) 

AT,— >oo 



which is equivalent to 



lim [(D c + m q )- x l5 + l5 (D c + m^ 1 } = , (27) 



m„— >0 



the chiral symmetry of the valence quark propagator (D c + m q )~ 1 in the mass- 
less limit. Therefore one can use this chirally symmetric (but nonlocal) D c to 
construct an exponentially-local lattice Dirac operator jU] 

D = D c (l + rD,)- 1 (28) 

which satisfies the Ginsparg- Wilson relation 

D l5 + l5 D = 2rD l5 D , (29) 

and its massive generalization [TT] 

D(m) = (D c + m)(l + rD c )~ l , (30) 

which is exactly in the same form of the effective 4D lattice Dirac operator 
(USD for any N s . 

Obviously, the optimal DWF action in the limit N s — > oo realizes the 
exact chiral symmetry (j2*7|) as well as the Ginsparg- Wilson relation ()29|). Fur- 
ther, for any finite N s , the action preserves the chiral symmetry optimally 
since the deviation of S op t from the sign function is the minimum, and has a 
theoretical upper bound [12], 



Ay = max 

v y^o 



ytQ2 y 

1 °opt 1 



<^^ = B z (N.,b) (31) 



where A is defined in (TL3|) . a function of iV s and b = ^max/^min- Thus, for 
any given gauge configuration, one can use the theoretical error bound to 
determine what values of N s and b (i.e., how many low-lying eigenmodes of 
should be projected out) are required to attain one's desired accuracy in 
preserving the exact chiral symmetry. In practice, with N s = 32, one should 
have no difficulties to achieve Az < 10~ 12 for any gauge configurations on a 
finite lattice (say, 16 3 x 32 at (3 = 6.0). The theoretical error bound can be 
estimated as 



B z (N s ,b) = A(b)e- c ® Ns 



where 



c(b) ~ 4.27(45) ln(6)-°' 746 ( 5 ) , 
A(b) ~ 8.11(1) r - 0091 ^ ln(6) 0042 ( 3 ) . 

Note that D ()28|) (N s — > oo) times a factor r = l/2mo is exactly equal to 
the overlap Dirac operator [T3*l I14j. 



D = - l+ 75 -^ . (32) 
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This implies that D is topologically-proper (i.e., with the correct index and 
axial anomaly), similar to the case of overlap Dirac operator. For any finite 
N s , r times D{m q = 0) (f2o"|) is exactly equal to the overlap Dirac operator D Q 
with (i/^) -1 / 2 approximated by Zolotarev optimal rational polynomial, (J1(J)) 
(N s = odd), or (N a = even). 

Athough the effective 4D lattice Dirac operator of the optimal DWF ac- 
tion is (^-invariant and preserves optimal chiral symmetry for any finite N 8 , 
the correlation function of quark fields q(x) and q(y) is not equal to the quark 
propagator (D c + m q )~^, nor the inverse of the effective 4D lattice Dirac oper- 
ator, D(m q )~y. This is unsatisfactory from the theoretical viewpoint. In this 
paper, I introduce two new boundary layers to the optimal DWF action, and 
derive the generating functional for n-point Green's function of quark fields 
q and q. Then the bilinear (and higher-order) quark operators manifest their 
continuum discrete symmetries, besides the as-invariance and optimal chiral 
symmetry for any finite N s . With 2D abelian background gauge fields, I com- 
pute the index, the axial anomaly, and the fermion determinant (ratio) of the 
effective lattice Dirac operator, for the optimal DWF and the conventional 
DWF respectively. All quantities of the optimal DWF are in good agreement 
with the exact solution, but the fermion determinant (ratio) of the conventional 
DWF disagrees with the exact solution, for topologically nontrivial sectors. 

The outline of this paper is as follows. In Section 2, we derive the generating 
functional for n-point Green's function of quarks fields in the optimal DWF 
as well as the conventional DWF, and discuss the issues pertainent to finite 
N s and a 5 . The discrete symmetry (V-A structure) of the left-handed quark 
current in effective weak Hamiltonian is shown explicitly. In Section 3, we 
compute the index, the axial anomaly, and the fermion determinant of the 
effective lattice Dirac operators in 2D abelian background gauge fields. In 
Section 4, we summarize and conclude with a discussion on the difference 
between the realization of exact chiral symmetry via the domain-wall fermions 
and the Ginsparg- Wilson lattice Dirac operator. In the Appendix, we derive a 
formula for computing the propagator of the effective 4D lattice Dirac operator 
(|23|) of the optimal DWF (0). 

2 Generating functional for n-point Green's 
function 

In quantum field theory, one of the fundamental quantities is the generating 
functional Z[J] (W[J) = \nZ[J]) for (connected) n-point Green's function. 
Once it is derived, one can obtain any (connected) n-point Green's function 
by successive differentiations with respect to the source J. Further, the gener- 
ating functional for one-particle irreducible vertices, r [</>], can also be obtained 
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through the Legendre transform T[(p] = W[J] — J dxJ(x)(j)(x). In this Section, 
we derive the generating functional Z\J] = e w ^ for n-point Green's function 
of quark fields q and q constructed from the boundary modes of domain-wall 
fermions. 

The generating functional for n-point Green's function of quark q and q is 
defined as 

p W[J,J\ = 7lT J ] = fJdU]\d#^ 

L ' J J[dU][d#ip)[d$\[d(j))e- A <>- A f-^ K ' 

where A g is the gauge field action, Af is the domain- wall fermion action, A p f 
is the corresponding pseudofermion action, and J and J are the Grassman 
sources of q and q respectively. 

We first derive Z[J, J] = e w ^ J '^ for the optimal DWF action, then show 
that it is invariant for any a$, and preserves the chiral symmetry optimally for 
any finite N s . Then we derive Z[J, J] for the usual DWF action, and discuss 
its problems due to finite iV s and 05. 



2.1 Optimal domain- wall fermions 

First, I generalize the optimal domain fermion action {I} by appending two 
boundary layers at s = and s — N s + 1, 

Ns+l 

A' f = ^2 ^2^(x,s){(uj s a 5 D w (x,x') + 5 X}X/ )8 StS/ 

s,s'=0 x,x' 

-{Sx,x> - uj s a 5 D w (x, x'))(P-6 s/ , s+ x + P + 5 s /, s _i)}^(a; / , s') 

N s +1 

= E J2^ s ) v '(^ s ^' s 'M x '^') (34) 

s,s'=0 x,x' 

with boundary conditions 

P+il>{x,-l) = -rm q P + ij(x,N s + l) , r = (35) 

P_ip(x, N s + 2) = -rm q P_^(x, 0) , (36) 

where the weights at the boundaries are zero, i.e. u = u>n s +i = 0, and other 
weights are given by (jZJ) as before. 

The quark fields coupling to physical hadrons are constructed from the left 
and right boundary modes 

q(x) = P_^(x,0) + P + ij(x,N s + l) , (37) 
q(x) = $(x, 0)P+ + $(x, N s + 1)P_ . (38) 

It will be shown that the fermion determinant of the new action (|34jl is 
exactly equal to that (fTHj) of the original action ((TJ). Thus the effective 4D 
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lattice Dirac operator is exactly equal to ([25)1. after the pseudofermion action 
(with m q = 2mo) is incorporated. The purpose of adding two boundary layers 
at s = and s = N s + 1 is for the construction of the quark fields (pT7 j) -(|38 |) 
such that the correlation function of q(x) and q(y) is equal to the valence quark 
propagator, i.e., r(q(x)q(y)) = (D c +m q )~^ in any background gauge field, and 
the generating functional Z[J, J] can be derived. 

Now we derive the generating functional Z [J, J] for the optimal DWF (|34|) , 
according to the formula (jH3J). 

Using 7sP± = ±P±, P+ + P- = 1, and H w = ^D w , we can rewrite (|3*4*|) as 

A' f = (rm q + 1)^075^+^+1 

-4> ol5 (P^ + P+^ Na+1 ) + ^075(^-^1 + P+^Po) 

Ns 

{^sJsi^sasHm - l)(P-1ps + P+lps-l) + 1psl^s^>H w + l)(P_Vs+l + P+^s)} 

s=l 

-^N 8 +ll5(P-^N s +l + P+^N S ) + ?/>7V s +l75(P-^0 + P+^N s +l) 

-(rm 5 + l)^ JVs+ i7 5 P^ (39) 

where all indices are suppressed except the flavor index in the fifth dimension. 
Next we define 

7? = P-ipo + P+ipN a +i , rj = -^o75 , (40) 
r] s = P-ip s + P + ip s -i , fj s = Tp s j 5 (u s a 5 H w - 1) ,s = 1,-- -,N S (41) 
Vn s +i = P-*Pn s +i + P+^N S , VN.+1 = -^w s +i75 , (42) 

and the inverse transform 

4>q = P-T] + P+% , vp = -r]o75 , (43) 
ip s = P-T) s + P+ris+i , ^ s = fj s (uj s a 5 H w - 1)^75 , s = 1, • ■ ■ , iV s (44) 
V'JVs+i = P-Vn s +i + P+% , ^jv s +i = -?7jv s +i75 • (45) 

Then the action (}39|) can be expressed in terms of {rj, fj} fields 
A' f = f}o(P--rm q P + )r} -f}orii + ^{f} s r} s -f) s T~ 1 r) s+ i} 

8=1 

+Vn s +iVn s +i ~ Vn s +i(P+ ~ rm q P_)r] (46) 

where 

T s = 1 = U ' a * Hw , a = l,...,JV.. 
1 + uj s a 5 H w 

Thus the fermionic integral in the numerator of (J33)) can be written as 

J [#] [rf^]e-' 4 /^^ 1+j9+9 " J = J J [dfj] [ dr] ] e -^v]+Jvo-v P + J+v Ns+1 P-J (47) 



where 

Jq + qJ = j(P-ipo + P + ip Nt+1 ) + $ Q P++<ii; Nt+1 PJ)J 

= Jr] -f] P + J + f] Ns+1 P_J (48) 

and J is the Jacobian of the transformation, 

N a +1 N a 

J= J] det [ l5 (u s a 5 H w - 1)] = [] det fr B (uj s a 5 H w - 1)] (49) 

S=0 8=1 

where u = ou Ng+ i = has been used. 

Now using the Grassman integral formula 

dxd X e- mxW * v = e* M ~ lv detM (50) 

one can easily evaluate the Grassman integrals in ([47)1 . by integrating (r) s ,fj s ) 
successively from s = N s + 1 to s = 0. 

Explicitly, after integrating (77^+1, f}N a +i), (|H|) becomes 

N. 



J / nt^*]^] ex P {~Vo(P- - rm q P + )rj + 7/0771 + Jr] - fj P + J 

N s -1 

~ (VsVs - fj s T-\ +1 ) - f] Ns r] Ns 
8=1 

+^ s T^[(P + - rm ? P_) % + P_ J] 
then integrating (t)n s ,Vn s ), the remaining integral becomes 

^ / II [ rf? 7^[^] ex P - rm q P + )r] + 77Q77! + Jt? - f] Q P + J 

J s=0 

N a -2 

- (VsVs - VsT^Vs+i) ~ *7jv.-i»7jv.-i 
8=1 

+^ s -iT^ 1 _ 1 T^ s 1 [(P + - rm q P-) Vo + P_ J]} 

Subsequent integrations over (r] Ns _ 2 ,f] Ns _2) up to (771,771) are similar to the 
above integration, and the result is 

jjm\im} exp {-MP. - r m ,P +)m + J„„ - ,„P+J 

+Vo II ^[(^H- - rm q P_) Vo + P_J] 



8=1 
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Finally, integrating (770, 770) gives 




(P_ - rm q P + ) - \{ T; l {P + - rm q P. 



s=l 



X 



-1 



p + + n T-'p-j \p. - n tt^+j ^ 

= K det[r(D c + m,)] exp { Jr~ l (D c + m g ) _1 j} 
where the simple identity 

— i . 

1 + 7 5 S opt 



rm n 



-1 



J 



-p+ + n T s ip - [ p --u t s i p+ 



s=l 



8=1 



1 - IhSopt 



rD r 



(51) 



(52) 



has been used, and K, S opt and rD c are defined in (j2DJ), (|2*2^1 and ({21} respec- 
tively. Note that any two transfer matrices commute with each other (T s T s i = 
T s fT s ), thus D c is invariant for swapping any two weights in {uii, ■ • • , ljn s }- 

Now setting J = J = in (|51|). we obtain the fermion determinant of the 
optimal DWF action 



detX>'(m,) = K det[r(P> c + m q )\ 



(53) 



where T>'{m q ) is the 5-dimensional lattice DWF operator defined in (|34jl . while 
rPv c is the 4-dimensional lattice Dirac operator defined in (|21|) . which becomes 
chirally symmetric in the limit iV, — > 00, the salient feature of domain- wall 
fermions. 

Obviously, if one repeats above steps to derive the fermion determinant of 
the original optimal DWF action (JTJ), then the result is exactly equal to (J53j) . 
The vital observation is that the transfer matrix T s at s = and s = N s + 1 
is equal to the identity (since uq = ^tv s +i = 0), thus 



N, 



N a +1 

n Ts^m 

S=0 8=1 



-1 



Explicitly, the transformation (j4"3|l - (|4"5j) is replaced by (fTtj jl - ffTSj) and the action 
(JTJ) can be expressed as 



N s -l 



A f [r], 77] = 7/1 (P_ - rrriqP^ - rj^ % + {naVs ~ VsT s V+i} + 

s=2 

Vn s Vn 8 - f]N 8 T N l(P + - rm q PJ)Tn . 
Thus the fermion determinant 



detV(m q 



[#][#] e^' 1 ^ = J 



-•Af[n,v] 
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can be evaluated by integrating {r] s ,f] s } successively from s — N s to s — 1, 
and the result is given in (j!9|) . 

As discussed in the Introduction, one cannot identify (D c + m q ) as the 
effective 4D lattice Dirac operator for the internal quark fields, since D c is 
nonlocal in the limit N s — > oo, and det(D c + m q ) does not reproduce 
det[7 A1 (9 M + iA IJL ) + m q ] in the continuum limit (a — > 0). The resolution is to in- 
troduce a set of pseudofermion fields {0, 0} similar to the domain-wall fermions 
but obeying Bose statistics. Explicitly, the action of the pseudofermion fields 
(with m q = 2mo) corresponding to the optimal DWF action (jH4j) is 

N a +1 

A' pf = ^2 ^2<j>{x, s){(uj s a 5 D w (x, x') + 5 XtX >)8 StS > 

s,s'=0 x,x' 

-{S x ,x> ~ uj s a 5 D w (x,x'))(P-5 s , iS+1 + P + 5 s / iS _ 1 )}0(x / , s') (54) 

with boundary conditions 

P+<f>{x,-l) = -P +( f>(x,N s + l) , 
P^(x, N s + 2) = -P-<p(x, 0) , 

where ujq = ujn s +i = 0, and other weights are given by (JTj) as before. Thus the 
integrations over the pseudofermion fields give 

J [d4>][d(j)]e' A ^ = n Ns+2 K~ l det(l + rD^ 1 (55) 

Substituting (|5*T|) and (|53j) into we obtain 

wu,J] _ J[dU]e- A °det[{D c + m q ){l + rDj- 1 ] exp {.Jr-\D C + m q )- l j) 
C ' ~ J[dU]e- A 9det[(D c + m q )(l + rD c )- 1 ] (56) 

This is one of the main results of this paper. 

Now it is clear that the most crucial quantity in (|5fi|) is D c which becomes 
chirally symmetric in the limit N s — > oo. At any finite N s , D c ()21|) preserves 
the chiral symmetry optimally since S opt ()22j) is equal to the optimal rational 
polynomial of the sign function Q. Moreover, u s a 5 (JJJ) is independent of a 5 , 
since it only depends on N s , \ min and A^ aa ./A^ in . It follows that T s , S opt , and 
D c are all independent of 05. Therefore the generating functional Z[J, J] for the 
optimal domain- wall fermion action (JMj) is a 5 -invariant, preserves the chiral 
symmetry optimally for any finite N s . Further, it also manifests the continuum 
discrete symmetries, since the quark propagator (D c + m g ) _1 (coupling to J 
and J) transforms like [j^(d p + iA^) + m q Y l . 

From (jnnj), it is clear that the effective 4D lattice Dirac operator for the 
internal quark fields (dressed with pseudofermions) is 

D(m q ) = (D c + m q )(l + rD.y 1 = m q + (m - m q /2)(l + j 5 H w R z (H w )) (57) 
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which is exactly equal to that (|25|) of the original action (fT|). 

Now any connected n-point Green's function can be obtained by differen- 
tiating W[J, J] with respect to J and J successively In particular, the quark 
propagator is 

5 2 W[J, J] 

r{q{x)q{y)} = - r- 



5J(x)5J(y) 



j=j=o 



J[dU}e^det[D(m q )} (D c + m q )-\x,y) 

J[dU]e-^det[D(m q )) 1 j 

which, in a background gauge field, becomes 

r{q{x)q{y)) = (D c + m q )~ l (x, y) , (59) 

and it goes to the Dirac propagator [j^(d lx +iA^)+m q )~ 1 (x, y) in the continuum 
limit. 

Now it is straightforward to generalize the optimal DWF action (|34|) to 
lattice QCD with rif quark flavors, with quark fields defined as 

q f (x) = V^[P-^f(x,0) + P + ^ f (x,N s + l)} , f = l,---,n f , (60) 
q f (x) = ^f(x,0)P + + ^ f (x,N s + l)P_] , (61) 

where the normalization factor y/r = l/\J2m has been included. Then the 
generating functional ()56|) can be rewritten as 

J[dC/]e-^n%idet[L>(m / )]ex P {E% 1 J/p c + m / )- 1 J / } 
Z[{J/ ' J/}] " ' J[dU]e-^U n f Udet[D(m f )} (62) 

where 

£)( m/ ) = [D c + m/)(l + rD c )~ l = m f + (1 - rm f )D 

Evidently, the generating functional (J62|) ensures that any bilinear (and 
higher-order) quark operator manifests the continuum discrete symmetries, 
since the quark fields ([BT} ]l -(|6T |l . both can be decomposed into left and right 
handed components with the usual chiral projectors, i.e. 

qj{x)=P-q f (x), qf(x) = P + q f (x), (63) 
qf(x) = q f {x)P+, qf(x) = q f {x)P_. (64) 

Thus the V-A structure of the left-handed quark current 

Ll^(z) = q fl (z) 1 ,(l- l5 )q f2 (z) (65) 
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in the effective weak Hamiltonian is preserved exactly. This is crucial for 
determining the indirect CP violation in the K —>■ tttt decay through the K°- 
K° mixing. The weak mixing matrix element 



(K»\s lfl (l - 75)^(1 - 7s)rf|^ U ) (66) 

(where s and d denote the strange and down quark fields respectively) can be 
extracted by evaluating the connected 3-point function (for all x and y and a 
range of 24 and 2/4) 

(J(x) 747 53(a;)L^(0)L^(0)%) 7475S ( 2 /)) connected 

^475777-7^77^7^(1 - 75)77777: x 



5J d (x) ^'°5J s (x)5J,(0) ,MV '°'SJ d (0) 

:7m(! - 7s) , j, n , , 77- x 7475 77777 Z[{J f , J/}] 



j=j=o 



6J S (Q) lo, SJ d (0) 5J d (y) '* lo 5J s ( y y 

— disconnected diagrams 

= \\ \dU\e~ A9 II det[D(m/)] tr[(D c + m d )^ xlil ^D c + m s )-^{l - 75)] x 

tr[(£> c + m 5 )-j7 A4 (l - 75 )(-D c + ™d) "j7475] , 
and the 2-point functions 

(rf(x) 7475 s(a;)s(0) 747 5rf(0)) 
= ^ / We-^ Il det [^K)] tr[(D c + m^^sPc + m,)j} W 5] , 

(S(0) 7475 C?(0)%) 747 5S(?/)) 

= 1 / W e ' A " U det l D ( m f)} tT l(D c + m s )-fa l5 (D c + m d )^ 7475 ] , 



where 



Z = J [dU]e- A » Y[det[D(m f )] 



f 

It is clear that there is no right-handed currents entering the weak mixing 
matrix element (|66|). and the SU(2)i gauge symmetry of electroweak inter- 
actions is preserved exactly. Had one not introduced the boundary layers at 
s = and s = N s + 1, then one might not have found any q and q such that Eq. 
floTIj) is satisfied and the usual chiral projectors (1 ± 75 )/2 can work for both q 
and q. Consequently its left-handed quark current may not manifest the V-A 
structure, and the SU{2) L gauge symmetry may be broken by the order 0(a) 
effect. Obviously, such left-handed quark current is not ideal for determining 
the K°-K° mixing. 
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In practice, the valence quark propagator (D c + nig)" 1 in a background 
gauge field can be computed via the 5-dimensional optimal DWF operator V 
()34jl . Using (jo^|) and the definition of quark fields (JH7j) - (JH%|) . one immediately 
obtains 



r 1 (D c + m q ) 1 (x,y) 
(q(x)q(y)} 



Z^l-^-^.O + -f+O s ,V s +lJ rr ; - ir - — (/+<V,0 + ^-<V,AT S +1, 



= £(P-5 s ,o + P+^+OP'- 1 ^, y, s')(P + <^,o + P-^,at s+ i) , (67) 

where s; y, s') can be obtained by solving the linear system 

V'(x,s;y,s')Y = 1 . (68) 

Alternatively, one can also obtain the valence quark propagator 

(D c + m,)- 1 = (1 - mg^fPrV?) " r] , (69) 

by computing D~ l (m q ) via the original optimal DWF operator (£[]), i.e., solving 
the linear system 

V(m q )Y = V(2m )B- 1 l , (70) 

where B and B" 1 are defined in (| 114)1 and ()115j) respectively. Then the solution 
Y gives 

D _1 (m ff )(a:,a:') = r(Sy)(a;,l;z',l) , (71) 

The derivation of ([70)) is given in the Appendix. 

Note that we can also compute the valence quark propagator (jHHJ) by solving 
the linear system 

m q + (mo - m q /2) (l + lb H w R z (H 2 w ))\ Y = 1 (72) 



D{m q )Y 

where D(m q ) is the effective 4D lattice Dirac operator (j57j). and Rz(H^) is 
the Zolotarev rational polynomial as defined in (j9~{)- (fTTj) . 

Since solving any one of above three linear systems (J72J), (fHHj) or (f7U|) yields 
the same quark propagator, in principle, it does not matter which linear system 
one chooses. However, in practice, one should pick the most efficient scheme 
for one's computational system (hardware and software). 
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2.2 Conventional domain-wall fermions 

Shamir's domain-wall fermion action j3] can be written as 

Afs = ^2^{x,s)[(a 5 D w (x,x') +S x y)S S: 

s,s'=l x,x' 



[P-5 S ^ S+1 + P + 5 s i )S ^i)5 X)X ,}il)(x\ s') 



= H ^^{x,s)V fs {x,s;x's')'4){x',s') (73) 

s,s'=l x,x' 

with boundary conditions 

P+V>0, 0) = -rm,? + ^(i, iV s ) , r = — (74) 
P_^(ar, N s + 1) = -rm q P^4>(x, 1) , (75) 

where Z)^ is defined in ()H), and m q is the bare quark mass. The quark fields 
are constructed from the boundary modes: 

q(x) = P^(x,l)+P+^(x,N s ) , (76) 
q(x) = $(x,l)P+ + $(x,N s )P-. . (77) 

The fermion determinant and the quark propagator (in gauge background) 
have been derived in Refs. PHJ-^Tj. Now we derive the generating functional 
for n-point Green's function of the quarks fields, according to the formula 
(|33|). The derivation is similar to that of the optimal DWF action in the last 
subsection. 

Using 75P± = ±P±, P+P- = 0, P + + P_ = 1, and H w = 7s-D„,, we can 
rewrite (f73|) as 

A fs = -{rm q + l)^i75(«5Pt«P- - l)P+^v s 

+^ l75 (a 5 P w P- - l)(P-^i + P+ip Ns ) + ^ 5 (a 5 H w P + + 1)(P_^2 + P+ipi) 

N s -1 

+E {*Psl5(a 5 H w P_ - 1)(P_V> S + P+ip s -i) + $ s j!>(a 5 H w P + + l)(P_^+i + P+^ s )} 

s=2 

+^PnM^H w P_ - l)(P_^7V s + P+^v s -i) + ^ N M^H W P + + 1)(P-^1 + P+lf>N.) 

-{rm q + l)$ Ns y 5 (a 5 H w P + + l)P-V>i (78) 

where all indices are suppressed except the flavor index in the fifth dimension. 
Next we define 

?7i = P_^i + P+ipN s , m = ^i75(«5Pt«P-- - 1) , 

r] s = P-ip s + P+ip s -i , fj s = ip s -f 5 (a 5 H w P^ - 1) ,s — 2,---,N s — l 

Vn s = P-4>N a + P+4>n s -i , Vn s = ^ Ns l5(a5HwP- - 1) , 
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and the inverse transform 

V>1 = P-7)! + P + 7]2 , $1 = %(a 5 P w P_ - 1) _1 7 5 , 

ip s = P_r? s + P + rj s+1 , ?/> s = f] s (a 5 H w P_ - l)~ lf y B , s = 2, • • • , N s - 1 
^jv 3 = P~Vn s + P+Vi , $n s = fj Ns (a 5 H w P_ - 1) -1 7 5 • 

Then the action (fTHj) can be expressed in terms of {77, f/} fields 

Ns-l 

A fs = ?7i(P- - rm q P + )r] 1 - fjiT' 1 ^ + ^ - VsT^rjs+i} 

s=2 



+Vn s Vn s - Vn s T \P+ - rm q P_y ni (79) 



where 



T = (1 + a 5 H w P + )-\l - a 5 H w P_) = \ ^ , (80) 

1 + O-5-n 

# = (2 + a 5 g w75 )- 1 g w = ^(2 + 7 5 a 5 ff w )- 1 = 75 . (81) 



At this point, it is instructive to compare the transfer matrix T (}80|) with 
its counterpart T s (|23|) in the optimal DWF. As it will become clear later, 
one of the difficulties in the conventional DWF is due to the fact that the 
denominator of H (|81|) is (2 + a^D w ) rather than a constant. 

Now the fermionic integral in the numerator of (|3*3*|) can be written as 

-A,,+Jq+qJ = T f [dfj}[dri}e- A ^ +Jm - f ' lP+J+ ^ T ~ lp - J (82) 



where 

Jq + qJ = J(P_ip i +P + ip Ns ) + (ip i P + +$ N P_)J 

= Jm + fii{^H w P_ - l)~ l P + J - f] Ns (a 5 H w P_ - iy l P_J 

= Jm-nP+j + t] Ns t- 1 p_j , 

and r is the Jacobian of the transformation, 

r = det [ lh {a h H w P„ - l)\ Ns . 

Again using the Grassman integral formula (|5()J). one can easily evaluate 
the Grassman integrals in (|82p. by integrating (j] s , fj s ) successively from s = N s 
to s — 1, similar to the steps in the last subsection. The final result is 



r det 



x 



(P_ - rm g P+) - T- Ns (P + - rm q P_ 
exp { J[(-P+ + T^P_)- 1 (P_ - T- Ns P + ) + rm q }- l j) 
K f det[r(pW + m g )] exp [jr-\D{ s] + m,)" 1 j} (83) 
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where 



K f = det[ l5 (a 5 H w P^ -l)] Ns det(- P + + T- N °P_) (84) 

rD is) = { _ p + T -N. p _ r i {p _ _ T -N sp+) = 1 + JsSm (85 ) 

1 - 756 (#) 

q( m (l + a 5 g)^-(l-a 5 g)^ 

1 J 1 + (1 + as/T)^ + (1 - a 5 H) N ' ' 1 J 

Again, since is nonlocal in the limit N s — > oo, one has to introduce 
pseudofermion fields (similar to the domain-wall fermions but obeying Bose 
statistics) with action (jTSj) (m 9 = 2m ) such that the effective 4D lattice 
Dirac operator for the internal fermion loops is exponentially-local for smooth 
gauge backgrounds, and it goes to [7 M (<9 M +iA fl ) + m q ] in the continuum limit 
a, 05 — ► 0. Then integrating the pseudofermion fields gives 

J [#][#] e- A "f* = tt n ° Kfdet(l + rD^Y 1 . (87) 

Substituting and (|87j) into we obtain 



e w[J,J] 



/[dE/]e-^det[(£>W + m ? )(l + rD^)- 1 ] exp { Jr~\D^ + m^ 1 j} 



/[dt/]e-^det[( J D£ ,) + m q ) (1 + r£> (sh 



which is exactly in the same form of (JSfiJl except replacing D c with Z)( s ) [i.e., 
replacing S opt (l22|) with S* (l86|) ]. However, such a difference has significant 
consequences, as we will demonstrate in the next section. 

In the limit N s — > 00, S(H) (jHEJ) becomes the sign function of H, 

lim S(H) 11 



N s ->oo yfjp 

thus ()85|) is chirally symmetric. However, it still depends on a 5 through 
H (JSTJ. In other words, the lattice spacing a 5 in the conventional DWF could 
be a relevant parameter 2 which must be taken to zero before its effective 4D 
lattice Dirac operator 

D (s) K) = (D^ + m q )(l + rD^y 1 (89) 

can behave properly. 

However, taking the limit 05 — > is rather tedious in practice. Therefore, 
it is vital to construct a domain-wall fermion action such that its generating 



2 As demonstrated by the tests in Section 3, the fermion determinant has a strong depen- 
dence on 05 for 2D topologically nontrivial gauge configurations. 
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functional Z[J, J] (thus any n-point Green's function) is independent of 05, for 
any a and iV s . 

This objective can be partially fullfilled if H (JHTj) is replaced by H w = 75-D™. 
Then, in the limit N s — > 00, S(H W ) (|5U|l becomes independent of a 5 , 

lim S(H W ) = -^L = Hw 

However, at finite N s , S(H W ) still depends on a 5 . 

Aiming at modifying the action (|75jl such that H is replaced by H w in the 
tranfer matrix T, one immediately observes that if all factors {a^H w P± ± 1) in 
the action (J7H|) are replaced by {a^H w ± 1), then T becomes 

r = . (90) 

This amounts to adding the following terms to the DWF action (J78j) . 

6Af s = -(rm q + l)ip 1 'y 5 a 5 H w P + ip Ns + 

+^j ll5 a 5 H w P + (P_^ 1 + P+i) Ns ) + i) llb ar H w P_(P_i) 2 + P+^i) 

N a -l 

+Y. {^sl^H w P + (P_i) s + P+^ s _i) + tfj s -f 5 a 5 H w P_(P_tjj s+1 + P+tps)} 

s=2 

+i>N s l<oa<oH w P + {P_i) Ns + P + 4> Na -i) + ip Na 'j 5 a 5 H w P_(P^tp 1 + P+^ Ns ) 
-(rm q + l^N^a^H^ipx , 

which can be simplified to 

6A fs = -rm q 'ip 1 j 5 a 5 H w P + ijj N3 + ^j^H^^ 

N s -1 

+ H {ipsl5a5H w P + ip s -x +il) s j 5 a 5 H w P-.ip s+ ij 

s=2 

+i>N a 'Y5a5H w P+if)Ns-i - rm (? ^ A r s 75a 5 if w P_^i . (91) 
Then adding (|9"T|) to (|75jl. one obtains 

^/s = XI H s){(a 5 -D w (x,x') + 8 XjX i)8 Sj8 > 

s,s'=l x,x' 

-(<**,*' - a 5 D w (x, x'))(P+5 s / iS _i + P_^, s +i)}^(x', s') , (92) 

which agrees with Borici's domain- wall action ^HJ . The effective 4D lattice 
Dirac operator of (J9*2~j) is the same as (|H9*j) except (the argument of D^) S(H) 
(|5U|) is replaced by S(H W ). Thus, at a 5 = 1, it (times a factor r = l/2m ) is 
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exactly equal to the overlap Dirac operator with (H^) l l 2 approximated by 
the polar approximation, since 



S(H V 



a ^ (i + W. Sr=i ^pfe) > N s = 2n+l (odd) 
a 5 H w I Er=r , N s = 2n (even) 



where 



sec 



N s \ 2 



d/ = tan 



N s V 2 



Obviously, for any finite N s , S(H W ) depends on a 5 . Also, it does not preserve 
chiral symmetry optimally, since S(H W ) is not the optimal rational approxi- 
mation of sga(H w ), as shown in Ref. [3]. 

A theoretical inconsistency of the Borici DWF action (|92|) is that the quark 
propagator {D' c +m q )~ l cannotbe expressed in terms of the correlation function 
of the quark fields constructed from the boundary modes (J75j) - (J77|) . i.e., 



r{q(x)q(y)} ^ (D' c + m q ) 



-l 

x,y i 



rD' 



1 + 1 5 S(H U 
1 - -y B S(H v 



in any background gauge field. Thus the generating functional (with localized 
quark sources J(x) and J(x)) for n-point Green's function of q and q is not 
well-defined, and cannot be derived. 

In this paper, a solution to this problem is presented in the last subsection. 
Explicitly, one can modify the action (|9*2*|) by adding two boundary layers with 
identity transfer matrix (i.e., without a 5 D w couplings) at s = and s — N s + 1, 
with boundary conditions ([35 p -(|36j h and defining the quark fields from the new 
boundary modes (|3*7 j) -(|3~5 }) . similar to the new optimal DWF (j31J). Then the 
generating functional for ra-point Green's function of q and q can be derived in 
the same way as that presented in the last subsection, and the result is in the 
same form of except replacing each T s (s = !,■••, N s ) with T' ()9(J|1. 



3 Tests 

In this section, we compute the index, the axial anomaly and the fermion 
determinant of the effective lattice Dirac operator of the optimal DWF ()34|) 3 . 
and the conventional DWF (|73*j) respectively, in two-dimensional background 
U(l) gauge fields. Our notations for the two-dimensional background gauge 
field are the same as those in Eqs. (7)-(ll) of Ref. [2*Uj . 

3 Note that both (JIJ and 1341) have the same effective lattice Dirac operator (|25p. as shown 
in Section 2. 
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The effective lattice Dirac operator of the optimal DWF can be computed 
from the 5- dimensional DWF operator ([34)1 . Setting m q = 2mo (rm q = 1) in 
floTj) . one has 



(l+rD c )-\x,y) 

= £(P-5 s ,o + P + 5 s , Ns+1 ){V'-\2m )}(x, s; y, s'){P + 5 s , fi + P-S^+i). 

s,s' 

Then the massless lattice Dirac operator can be obtained by a subtraction 

D = D c {l + rD c y 1 = \-{l+rD c y 1 , (93) 
and the massive one by 

D{m q ) = (D c + m g )(l + rD^ 1 = 1 + (m q - 1)(1 + r^)- 1 . (94) 

Similarly, we can obtain the effective lattice Dirac operator of the conven- 
tional DWF ()73j). with the boundary modes at s = 1 and iV s . 

In the following tests, we fix m = 1 (r = 1/2) in D w (£Q), and m q = 
(except for the fermion determinants in Table 2 with m q = 0.05). With N s in 
the range 10 < N s < 32, we first check that the eigenvalues of the effective 
lattice Dirac operators fall on a circle of radius one (l/2r = 1) centered at one. 
For nontrivial gauge backgrounds, the real eigenmodes (zero and two) have 
definite (+1 or —1) chirality and satisfy the chirality sum rule [20] 

n+ + N + = n_ + AL , (95) 

where n±(N±) stand for the numbers of zero (+2) eigenmodes with ±1 chiral- 
ity respectively. These real eigenmodes are topological, i.e., robust under local 
fluctuations of the gauge fields. For any (trivial or nontrivial) gauge back- 
ground, the complex eigenmodes must come in conjugate pairs, each with zero 
chirality. The Atiyah-Singer index theorem (n + — n_ = Q) holds for any gauge 
configuration (with topological charge Q) satisfying the topological bound [22] 

\f x \p(x)\ < 0.28 (96) 

where p(x) is the topological charge (density) inside a unit square (of area a 2 ) 
centered at x, 

I rxi+a/2 rx 2 +a/2 \ 

Pix) = - dy x dy 2 —F l2 {y) . (97) 

GT Jxi-a/2 Jx 2 -a/2 Z7T 

Thus, we conclude that the effective lattice Dirac operators of the optimal 
DWF and the conventional DWF both are topologically-proper for sufficiently 
large N s . 
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3.1 Axial anomaly 



Since the effective lattice Dirac operator D of domain-wall fermions (N s — > oo) 
satisfies the Ginsparg- Wilson relation (J29j) . its axial anomaly (for any N s ) can 
be defined as jTH] 

A L (x) = tr[ 75 (l - rD(x, x))\ = --tr[5(x, x)} , (98) 

where the trace runs over Dirac and color indices. Here the operator S becomes 
S pt (122) for the optimal DWF, and S §^ for the conventional DWF. 

For sufficiently large N s and for gauge configurations satisfying the topo- 
logical bound ()96|) . I checked that the sum of the axial anomaly over all sites 
is equal to the index of D and the topological charge Q 

J2A L (x) = n + -n_ = Q , (99) 

X 

for both the optimal DWF and the conventional DWF. However, this does not 
necessarily imply that Al(x) is equal to the topological charge density at each 
site. 

In the following, I compute the axial anomaly Al{x) for the overlap Dirac 
operator, the optimal DWF, and the conventional DWF respectively, and com- 
pare them with the topological charge density p(x) (fHTj) of the gauge back- 
ground on the torus. The deviation of the axial anomaly can be measured in 
terms of 

^=^El^(x)-a 2 p(x)| (100) 

where p(x) is the topological charge density inside a unit square (of 
centered at x. 

In a nontrivial U(l) gauge background with parameters Q — 1, hi — 0.1, 
h 2 = 0.2, Af ] = 0.3, A { 2 ] = 0.4 and n x = n 2 = 1 (as defined in Eqs. (7) and 
(8) in Ref. [20]) on the 8x8 lattice, the axial anomaly Al{x) and its deviation 
5 are computed for the overlap, the optimal DWF and the conventional DWF 
respectively. Here the 2D overlap Dirac operator is computed exactly via 
diagonalization. 

The results are : 

!4.75 x 10 -3 , overlap , 
4.75 x 10~ 3 , optimal DWF (N s = 16) , (101) 
5.44 x lO" 3 , conventional DWF (N s = 16, 32) . 

The slightly larger deviation of the axial anomaly of the conventional DWF 
indicates that its effective lattice Dirac operator is less localized than that of 
the optimal DWF. Note that there is no improvements by increasing iV s = 16 
to N s = 32. 

For all background gauge configurations we have tested, the optimal DWF 
always gives anomaly deviation (5) smaller than that of the conventional DWF. 
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3.2 Fermion determinant 



In quantum field theory, the fermion determinant det(D) is proportional to the 
exponentiation of the one-loop effective action which is the summation of any 
number of external sources interacting with one internal fermion loop. Thus 
it is one of the most crucial quantities to be examined in any lattice fermion 
formulations. Even though a (massless) lattice Dirac operator is topologically- 
proper (i.e., with the correct index and axial anomaly), it does not necessarily 
imply that its non-zero complex eigenvalues behave properly. This can be 
tested by computing its fermion determinant 



where det'(D) is equal to the product of all non-zero eigenvalues. Since the 
effective lattice Dirac operators of the optimal DWF and the conventional 
DWF are 75-hermitian, the eigenvalues of D are either real or come in complex 
conjugate pairs, and det'(D) must be real and positive. For Q = 0, we have 
n + + n_ = and det(D) = det'(D). For Q 7^ 0, we have n + + n_ 7^ 
and det(-D) = 0, but det'(D) still provides important information about the 
spectrum. In continuum, exact solution of fermion determinant in the general 
background U(l) gauge fields on a torus (L\ x L 2 ) was obtained in Ref. [2~T] . 

With nontrivial gauge backgrounds, we compute det'(D) for the overlap 
Dirac operator, the optimal DWF, and the conventional DWF respectively, 
and compare them with the exact solution in continuum. For simplicity, we 
turn off the harmonic part (hi — h 2 — 0) and the local sinusoidal fluctuations 
(A^ = A 2 = 0), and examine the change of det'(D) with respect to the 
topological charge Q. For these gauge configurations, the exact solution [2*T] 
is 



such that det'[D(l)] = 1. 

In Table 1, the fermion determinants det'(D) on the 8x8 lattice are listed. 
It is clear that the fermion determinant of the optimal DWF (N s = 10) highly 
agrees with that of the (exact) overlap Dirac operator, and both of them are in 
good agreement with the exact solution on the 8x8 torus. On the other hand, 
the fermion determinant of the conventional DWF (N s = 10, = 1) does not 
agree with the exact solution. Further, the discrepancies do not decrease as 
N s is increased. This suggests that the error is not due to finite N s , but the 



det(D) = (1 + e i7r ) (n++n - ) det / ( J D) 



(102) 




(103) 



where the normalization constant M is fixed by 
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exact 


overlap 


optimal DWF 


conventional DWF 


Q 






N 8 = 10 


N s = 10 
a 5 = 1 


N s = 32 
a 5 = 1 


iV s = 20 
a 5 = 0.5 


1 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


2 


2.82843 


2.79094 


2.79094 


1.44349 


1.44558 


2.08267 


3 


6.15840 


5.96482 


5.96482 


1.63969 


1.64629 


3.32119 


4 


11.3137 


10.7257 


10.7257 


1.61127 


1.62411 


4.46412 


5 


18.3179 


16.9800 


16.9800 


1.43207 


1.45172 


5.29955 


6 


26.8177 


24.3215 


24.3215 


1.18079 


1.20677 


5.71414 


7 


36.1083 


32.0978 


32.0978 


0.91741 


0.94853 


5.69848 


8 


45.2548 


40.5552 


40.5551 


0.67521 


0.71130 


5.37711 


9 


53.2732 


45.9271 


45.9269 


0.47955 


0.51653 


4.68732 


10 


59.3164 


50.6766 


50.6761 


0.32497 


0.36279 


3.91209 



Table 1: The fermion determinant versus the topological charge Q. The nor- 
malization constant is chosen such that det'D = 1 for Q = 1. The results on 
the 8x8 lattice are listed for the overlap Dirac operator, the effective lattice 
Dirac operators of the optimal DWF and the conventional DWF respectively. 
The exact solutions on the 8x8 torus are computed according to Eq. (|103|) . 

lattice spacing which enters the effective lattice Dirac operator (J%9"|) through 
S(H) (j8Hj) . If this is the case, then decreasing a$ would improve the results. 
This is demonstrated in the last column with = 0.5 and N s = 20, which 
gives results much better than those of the column with N s = 32, 05 = 1. 

In Table 2, the fermion determinants for massive fermion with m q = 0.05 
are listed, for the same nontrivial gauge backgrounds in Table 1. The massive 
overlap Dirac operator is [m q + (1 — rm q )D ], and the effective lattice Dirac 
operators of the optimal DWF and the conventional DWF are computed via 
(I94j) . In each case (column), the normalization constant is fixed such that 
det[D(m q )] — 1 for Q — 0. Clearly, the fermion determinant of the optimal 
DWF (N s = 10) is in excellent agreement with that of the (exact) overlap, 
while the conventional DWF disagrees with the overlap and the optimal DWF. 
Further, the discrepancies do not decrease as N s is increased. However, if 05 is 
reduced from 1 to 0.5, then the results are significantly improved, as shown in 
the last column. This demonstrates that the massive fermion determinant of 
the conventional DWF is also sensitive to a 5 , for topologically nontrivial gauge 
configurations. 

On the other hand, for 2D topologically trivial gauge configurations, the 
fermion determinant of the conventional DWF only has a weak dependence on 
a 5 . This can be demonstrated as follows. For simplicity, we turn off the local 
sinusoidal fluctuations (A^ = = 0), and examine the change of det(.D) 
versus the harmonic parts hi = h-2. In Table 3, the fermion determinants 
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overlap 


optimal DWF 


conventional DWF 


Q 




N 8 = 10 


N s = 10 
a 5 = 1 


N s = 32 
a 5 = 1 


N s = 20 
a 5 = 0.5 





1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1 


5.27834 


5.27834 


10.4991 


10.4943 


7.00364 


2 


0.73042 


0.73042 


0.75612 


0.75684 


0.72616 


3 


7.768 x 10~ 2 


7.768 x 10" 2 


4.289 x 10~ 2 


4.304 x 10~ 2 


5.782 x 10- 2 


4 


6.961 x 10~ 3 


6.961 x 10" 3 


2.106 x 10" 3 


2.121 x 10" 3 


3.893 x 10~ 3 


5 


5.497 x 10~ 4 


5.497 x 10" 4 


9.371 x 10" 5 


9.474 x 10" 5 


2.332 x 10~ 4 


6 


3.929 x 10~ 5 


3.929 x 10" 5 


3.883 x 10" 6 


3.935 x 10~ e 


1.289 x 10~ 5 


7 


2.588 x 10~ 6 


2.588 x 10" 6 


1.535 x 10~ 7 


1.535 x 10" 7 


6.830 x 10~ 7 


8 


1.633 x 10~ 7 


1.633 x 10~ 7 


5.621 x 10" 9 


5.794 x 10" 9 


3.407 x 10~ 8 


9 


9.235 x 10~ 9 


9.235 x 10" 9 


2.405 x 10" 10 


2.103 x 10" 1U 


2.156 x 10~ 9 



Table 2: The fermion determinant det[D(m q )] versus the topological charge 
Q for effective lattice Dirac operator with bare mass m q = 0.05 on the 8x8 
lattice. The normalization constant is chosen such that det D(m q ) = 1 for 
Q = 0. The results are listed for the overlap, the optimal DWF (N s = 10), 
and the conventional DWF with (N s = 10, = 1), (N s = 32, a*, = 1) and 
(N s = 20, 05 = 0.5) respectively. 

versus the harmonic parts (hi = h 2 ) on a 8 x 8 torus (lattice) are listed for the 
exact solution [21], the overlap, the optimal DWF, and the conventional DWF 
respectively. Clearly, all results are in very good agreement with each other. 
Further, the fermion determinant of the conventional DWF is insensitive to a^, 
as comparing the corresponding entries of the last two columns (N s = 10, a 5 = 
1) and (iV s = 20, a 5 = 0.5). 

Even though above tests are performed with 2D abelian gauge backgrounds, 
one expects that the a 5 effects might also occur in the quark propagator (as 
well as the fermion determinant), in lattice QCD with the conventional DWF. 
We are now computing the quark propagator for several different values of 
05, on the 8 3 x 24 lattice, for 100 gauge configurations. Then we will extract 
the pion mass as well as other observables, and study their dependence on 05. 
These results will be presented elsewhere. 

4 Concluding Remarks 

To formulate lattice QCD with domain-wall fermions, it seems to be inevitable 
to introduce parameters N s and a 5 (in the fifth dimension) to the theory. 
Presumably, only in the limit iV s — > 00 and a 5 — > 0, the effective 4D theory 
agrees with the target theory (QCD), up to the finite volume and discretization 
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exact 


overlap 


optimal DWF 


conventional DWF 


hi = h 2 






N s = 10 


N s = 10 

a 5 — 1 


N s = 20 
a 5 = 0.5 


0.0 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.1 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


0.2 


3.34414 


3.40518 


3.40518 


3.39253 


3.40883 


0.3 


5.71957 


5.97173 


5.97137 


5.92350 


5.98565 


0.4 


7.24506 


7.73701 


7.73701 


7.64986 


7.76372 


0.5 


7.73926 


8.34167 


8.34166 


8.23777 


8.37383 


0.6 


7.24506 


7.73701 


7.73701 


7.64986 


7.76372 


0.7 


5.71957 


5.97173 


5.97137 


5.92350 


5.98565 


0.8 


3.34414 


3.40518 


3.40518 


3.39253 


3.40883 


0.9 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.0 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 



Table 3: The fermion determinant versus the harmonic parts (hi = h 2 ), on the 
8x8 lattice. In each column, the normalization constant is chosen such that 
det(-D) = 1 for hi — hi — 0.1. Note that the fermion determinants satisfy the 
symmetries of the exact solution: (i) periodicities in hi and h 2 with period 1; 
(ii) reflection symmetries with respect to hi — 1/2 and h 2 = 1/2. 



errors. The relevant question is how to minimize the dependence on these two 
extra parameters. More importantly, if the dependence on 05 is too strong, 
then at nonzero 05, the effective 4D theory might differ significantly from the 
target theory, even in the limit N s — > 00. Therefore, it is desirable to eliminate 
a 5 completely and to preserve the chiral symmetry optimally for any N s . 

In the limit N s — > 00 followed by 05 — > 0, the fermion determinant (includ- 
ing the contribution of pseudofermions) of the conventional DWF gives the 
effective 4D lattice Dirac operator 



where 



D(m g ) = {D c + m q ){\ + rD^ 1 =m q + (l- rm q )D 



1-75^W)- 1/2 
D = Dcil + rDj-^moll + ^H^Hly 1 / 2 ], 

while the quark propagator in background gauge field is 

{q(x)q(y)) = (D c + m q )£ y . (104) 

Since the optimal rational approximation for sgn(H w ) = H^H"^)' 1 / 2 is 
H W R Z (H 2 ), where Rz{x) is the Zolotarev optimal rational polynomial for 
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x -1 / 2 (i.e., R^~ l ' n \x) for N s = 2n, or R% (x) for N s — 2n+ 1), the relevant 
question is whether one can construct a DWF action such that its effective 
lattice Dirac operator for finite N s is 

D(m q ) =m q + (m - m q /2)[l + ^H W R Z (H 2 W )\ , 

in view of that the S operator in the conventional DWF is a rational polynomial 
approximation of sgn(H). The solution to this problem is the original optimal 
DWF action [5 . The salient feature of the optimal DWF is that its weights 
{uj s } © not only give S opt , the optimal rational approximation of sgn(Hu,), 
but also eliminate 05 completely. 

However, the correlation function of q(x) and q(y) is not equal to the quark 
propagator (D c + m q )~^, nor the inverse of the effective 4D lattice Dirac op- 
erator, D[m q )~}y. In other words, either (D c + m q )~ x or D(m q )~ x cannot be 
expressed in terms of the correlation function of local quark fields q(x) and q(y), 
but q(x) and some "nonlocal" field variable Q(y) which depends on q(x) at all 
x. This posed a rather serious problem to the optimal DWF, since it implies 
that the basic concept of local fields q(x) and q(x) in quantum field theory has 
to be abandoned in the optimal DWF. Even if Q(y) might be exponentially- 
local, it is still unsatisfactory in comparison with the exactly local quarks fields 
q and q in the conventional DWF. 

The solution to this problem is presented in this paper: Appending two 
boundary layers at s = and s = N s + 1 with u = wat s +i — 0. Then the 
correlation function of (local) quark fields q{x) and q(y) constructed from these 
new boundary modes (pT7 j) -(p?8 j) satisfies (jlU4|) . Consequently, one can derive 
the generating functional for n-point Green's function of quark fields (with 
local Grassman sources J and J) 

W[J,J] _ J[dU}e- A *det[(D c + m q )(l + rD c )- 1 } exp { J(D C + m^j) 
C ' ~ /[d£7]e-^det[(D c + m q )(l + rD c )- 1 } ' 

Then the bilinear (and higher order) quark operators manifest their contin- 
uum discrete symmetries at finite lattice spacing, since the chiral projectors for 
both q and q are the usual P± — (1 ± 7s)/2. In particular, the V-A structure 
of the left-handed quark current in the effective weak Hamiltonian is preserved 
exactly at any finite lattice spacing, which is crucial for investigating the weak 
matrix elements in lattice QCD, e.g., the CP violation in the K — ► tttt decay 
through K°-K° mixing. 

From the 2D tests in Section 3, we see that the fermion determinant (ra- 
tio) of the effective lattice Dirac operator D of the conventional DWF does 
not agree with the exact solution for any topologically nontrivial background 
gauge field, even though D can capture the topology of the gauge background 
(i.e, reproducing topological zero and nonzero real eigenmodes with definite 
chirality, satisfying the Atiyah-Singer index theorem on a finite lattice, and its 
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axial anomaly in agreement with the topological charge density) . It is interest- 
ing to see that a lattice Dirac operator satisfies the Ginsparg- Wilson relation, 
75-hermiticity, doubler-free, exponential-locality, and the correct continuum 
(free fermion) behavior, but could not reproduce the proper fermion determi- 
nant for topologically nontrivial gauge backgrounds, due to its dependence on 
05(7^ 0). This warrants further studies on effects in lattice QCD with the 
conventional DWF (e.g., to investigate to what extent the pion mass or other 
observables depending on a 5 ). 

Finally, it is instructive to compare the quark fields in lattice QCD with 
Ginsparg- Wilson Dirac operator to the exactly local quark fields in the (op- 
timal/conventional) DWF. In general, the quark action of lattice QCD with 
Ginsparg- Wilson Dirac operator reads 

A ^J2^( x ) D ^y( m i)i(y) ' ( 105 ) 

x,y 

where 

D(m q ) = {D c + m g )(l + rD.)' 1 = m q + (1 - rm q )D, D = D c (l + rD c )~ l , 

and D satisfies the Ginsparg- Wilson relation (J2HJ). Here we assume that D 
is doubler-free, exponentially-local for sufficiently smooth background gauge 
fields, and has the correct continuum (free fermion) behavior. Then the quark 
propagator in background gauge field is 

<9(*MV)> = f[dq][dq]e-A> = ^- K) ' (1 ° 6) 

But this cannot be the valence quark propagator since it breaks the chiral 
symmetry even in the massless (m q = 0) limit. Thus the quark propagator has 
to be extracted as 

(D c + m q )zl = (1 - rm q y l [D-l{m q ) - r5 x , y ] . 

In other words, one has to replace either q(x) or q(y) in the quark bilinear of 
(jl06J) with one of the followings: 

q(x) -> q(x) = £(1 + rD c )-] z q{z) = £(1 - rD) X:Z q(z) , 

z z 

or 

q(y) - f(2/) = E^)( 1 + ^)^ = E^)( 1 -^ky 

z z 

Note that both above transformations are exponentially-local. Suppose one 
replaces q(x) with q(x), then the quark propagator (in background gauge field) 
becomes 

(q(x)q(y)) = (D c + m q )^ y . (107) 
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Even though the concept of local field is only violated exponentially by q(x), 
it is still unsatisfactory from the theoretical viewpoint, in comparison with the 
exactly local q and q in the (optimal/conventional) DWF. To impose the exact 
locality of q, one can rewrite the quark action as 



which, however, is also unsatisfactory since D c is nonlocal. In general, it seems 
to be impossible to formulate lattice QCD with any Ginsparg- Wilson lattice 
Dirac operator such that the quark fields and the action both are consistent. 
However, if the Ginsparg- Wilson lattice Dirac operator (e.g., the overlap Dirac 
operator with Zolotarev rational approximation) is equal to the effective 4D 
lattice Dirac operator of certain DWF action (e.g., the optimal DWF), then 
one can resort to the DWF action, in which the quark fields q and q (e.g., ([37)1- 
(J2H1)) are well-defined in terms of the boundary modes, with exact locality as 
well as the usual chiral projectors P± — (1 ± 7s)/2. 

Appendix 

In this appendix, we derive a formula for computing the inverse of the effective 
4D lattice Dirac operator, with the original optimal DWF operator (£[]). 
From (|T^|) . we have 

det[P(m g )P" 1 (2m )] = det[r(D c + m q )(l + rD c y l ] = det[rD(m q )) (109) 
or equivalently, 



where T> is the optimal DWF operator defined in while D(m q ) is the 
effective 4D lattice Dirac operator. It follows that there may exist an operator 
Q such that 



{QV- 1 (m q )V(2m )Q- l }(x,l;x',l) = r^D-^m^ix,^) . (Ill) 



Now we want to find an explicit expression of Q. Using the Grassman 
integral formula 



{QV- 1 (m q )V(2m )Q- 1 }(x, 1; x', 1) 

A * s / [#][#KW(*, l){^V(2m )Q- l }(x', l)e' A f . (112) 
detZJ (rrin) J 



A = Yl <l( x )( D c + m q ) x , y q(y) , 



(108) 



det[V-\m q )V(2m )} = det^ 1 D' 1 (m q )] 



(110) 




we rewrite the l.h.s. of (jlllj) as 
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From the transformation ([16 p -(|18j h one immediately sees that the optimal 
DWF action (0) can be rewritten as 



Af = Y^s{x)V{m q ){x,s;x's')^ s >{x') 

s,s'=l x,x' 

= E I! fj s (x){(tta 5 H w - l)~ 1 -f 5 V(m q )B~ 1 }(x, s; x', s')r) s ,(x') 



s,s'=l x,x' 

Ns-l 

Vi(P- — rm qP+) r li ~ ViT{ l i]2 + {VsVs ~ VsT~\ + l) 



s=2 



where 



and 



+Vn s Vn s ~ Vn s T Ns (P+ - rm 9 P_)^ (113) 



Q(x, s; x', s') = u s 5 Si s'5 X:X > 



B(x,s;x',s') = 5 X}X >(P-5 SiS , + P+5 S _ M , (114) 
B~ 1 (x,s;x',s') = 5 x>x ,(P-5 SjS , + P+5 s+ljS ,) , (115) 

with periodic boundary conditions in the fifth dimension. 

Now it is clear that Q = B would satisfy (jlllj) . This can be seen as follows. 
Using the transformation (|16j) - (jl8j) . the r.h.s. of (|112j) can be written as 



J 



detV(m q 



[df\}[drj\e- A ^ 1 {x){fi{Q.a b H w - l)- 1 l5 V{2m )B- 1 }(x', 1) (116) 



where J is the Jacobian defined in (J4*9~j) . Using ()113p . the last expression in 
the integrand of (jllfij) can be simplified as 

{fj{na 5 H w - l)- 1 l5 V(2m )B- 1 }(x', 1) 
= f}x(x')(P- - P+) - {fj Ns T£}(x')(P + - P-) . (117) 

Thus ()116|) becomes 

3 frr-in l -J, f / \ - / /\ / \ r - m-li f l\ 



, , , l-wjl-wj- -r]x{x)r]x{x') - Vi(x){VN a T Ns }(x') 75 , (118 

detP(m 9 ) J 1 J 

which can be integrated successively from (r] Ns ,f] Ns ) to (771,771). Explicitly 

[dff\ [dr l ]e' A f'q 1 {x){f] N T N l s }(x') 75 
=/[^ 1 ][d ?7l ] e -^[^-- rm « F ^-nf4^ 1 (^-™^-)]m rnix^fh IIt^x^ts , 
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which is substituted into ()118|) to give 
J 



detV(m q 



y - [cZ77 1 ][^ 1 ]e-^[( p --^ p +)-nf= 1 ^ 1 (^-^-)] 



in 



s=l 



{.*) 75 



(P_ - rm 9 P + ) - J] " rm q PJ y 



s=l 



= r- 1 [(D c + m q )- 1 (l + rD c )](x 1 a/) = r^D^im^x, x') 

where ()49|h (fTTJJ) and (Jo"2"|) have been used. 
Therefore, we have 

{BV- l (m q )V(2m )B- 1 }(x, 1; a/, 1) = r" 1 £) _1 (m,)(x, x') 

and D~ l (m q ) can be obtained by solving the linear system 

V(m q )Y = V(2m )B- 1 l . 

The solution of (|12(J|) gives 

i)- 1 (?n J )(z,i')=r(By)(i l l;i , ) l) , 

and the valence quark propagator (J69j) . 



(119) 



(120) 



;m) 
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